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ABSTRACT 

Several indirect evidences indicate a magnetic origin for the non-thermal 
width of spectral lines observed toward molecular clouds. In this letter, I suggest 
that the origin of the non-thermal width of carbon recombination lines (CRLs) 
observed from photo-dissociation regions (PDRs) near ultra-compact H II regions 
is magnetic and that the magnitude of the line width is an estimate of the 
Alfuen speed. The magnetic field strengths estimated based on this suggestion 
compare well with those measured toward molecular clouds with densities similar 
to PDR densities. I conclude that multi-frequency CRL observations have the 
potential to form a new tool to determine the field strength near star forming 
regions. 

Subject headings: ISM:magnetic fields - radio lines:ISM - ISM:lines and bands 
- ISM:molecules - MHD - H II regions 



1. Introduction 



Observations of high-density molecular line tracers toward ultra-compact Hn regions 
(UCHs) reveal that these Hn regions are embedded in dense (^10 5 cm -3 ) molecular clouds 
(eg. Churchwell, Walmsley & Cesaroni 1990; Churchwell 2002 and references there in). Far 
ultra-violet (FUV; 6 to 13.6 eV) radiation from massive stars within the UCHs heats the 
dense molecular material in this interface, producing a Photo Dissociation Region (PDR) 
(see Hollenbach Sz Tielens 1997). Carbon recombination lines (CRLs) from such regions have 
been detected toward a large number of UCHs, establishing the presence of dense PDRs near 
most UCHs (Roshi et al. 2005a). 
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The observed width of the CRLs from PDRs associated with UCHs is typically between 
4 and 8 km s _1 . The gas temperature in the PDR, obtained by non-LTE modeling of CRL 
emission at multiple frequencies, is in the range 200 to 1000 K (Roshi et al. 2005b; Garay 
et al. 1998; Natta, Walmsley & Tielens 1994). Thus the expected thermal contribution to 
the line width is at least a factor of two smaller than the observed width (see Fig. [1]). This 
larger observed width indicates that the CRL width is dominated by supersonic motions. 
In fact, a similar situation exists in molecular clouds. It has long been recognized that the 
observed width of spectral lines from molecular clouds have a non-thermal component and 
this component is often supersonic (eg. Barrett, Meeks and Weinreb 1964). The origin of the 
non-thermal component of the line width is now attributed to either "turbulence" (Morris 
et al. 1974; Zuckerman & Evans 1974) or Alfuen waves (Arons & Max 1975; Mouschovias 
1975; see also Shu, Adams & Lizano 1987 and references there in). 

Detailed theoretical models for both "turbulence" and Alfuen waves in molecular clouds 
are yet to be worked out. But "turbulence" without any magnetic field is a less plausible 
proposition since observations demand supersonic turbulence and the rapid dissipation in 
shocks cannot sustain such turbulent motions to time scales larger than the free-fall time 
(eg. Shu et al. 1987). The decay of Alfuen waves are slow compared to supersonic turbulence 
though simulations show that they may also need a driving mechanism (MacLow 2003). Thus 
Alfyen waves appear to be a more tenable explanation for the observed line widths (Arons 
& Max 1975; Mouschovias 1975). 

Observational data also seem to indirectly support the magnetic origin of the line width. 
Analysis of the width of spectral lines observed toward molecular clouds shows power-law 
relationships between the non-thermal component of the velocity dispersion, the cloud size, 
density and the magnetic field strength (Larson 1981; Troland & Heiles 1986; Myers & 
Goodman 1988a, 1988b; Mouschovias & Psaltis 1995; Crutcher 1999). These relationships 
are expected in self-gravitating, magnetically supported clouds since the Alfuen waves in 
such clouds affect the observed line width (Mouschovias 1987; Myers & Goodman 1988a, 
1998b). 

In this letter, I use the observed non-thermal widths of CRLs to estimate the magnetic 
field strength in PDRs near UCHs. The field strength is estimated by assuming that the 
non-thermal velocity dispersion of carbon lines give an estimate of the AIKen speed in the 
PDR. This assumption is based on a comparative study of the characteristics exhibited by 
the CRL and molecular line data. The data and the comparative study are presented in 
Sections [2] and [3] respectively. The magnetic field strength is obtained by combining the 
velocity dispersion of CRL and the density derived from modeling carbon line emission (see 
Section H]). To our knowledge the only attempt to compute the magnetic field using CRLs 
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was made by Vallee (1989), where he had applied a model for shocks at the edge of molecular 
clouds to deduce the field strength. 

2. Carbon recombination line data 

For the present analysis I use CRL data toward UCHs since the PDR thickness is small 
in these cases and therefore large scale motions (such as outflows) may not affect the line 
width (see Section[5]). Table. [1] summarizes the CRL data and results obtained from modeling 
this data toward 14 UCHs. Listed are the source name, observed CRL transition, FWHM of 
CRLd3, PDR gas temperature, carbon ion density, number density of hydrogen atoms in the 
PDR, line-of-sight extent of the PDR (Ly) and the estimated magnetic field strength (see 
Section SD. 

The quality of the CRL data obtained toward all 14 sources is good; the uncertainties are 
mostly in the parameters derived by modeling line emission. The frequencies of the selected 
CRLs are below ~ 15 GHz. As inferred from modeling, carbon line emission is dominated by 
stimulated emission at these frequencies ( Natta et al. 1994; Roshi et al. 2005a). The data 
toward the first 7 sources in Table [1] were obtained using interferometric observations. The 
line widths for the first six sources were obtained from profiles averaged over regions where 
lines were detected. The properties of the PDR toward W3A were not well constrained; 
I have taken representative model parameters from Kantharia, Anantharamaiah & Goss 
(1998). Data toward the remaining 7 sources were obtained with the Arecibo telescope. 
Modeling of CRL emission toward these sources was limited by: (a)poor angular resolution, 
(b)flux density calibration error and (c) detection of line only near 9 GHz. Thus modeling 
this data set has given only constraints on the physical properties of the PDR. The PDR 
gas temperature of 500 K for the Arecibo sources and W48A, W49G and W49J listed in 
Table [1] is a representative value. The carbon ion density and L\\ obtained from modeling 
depend non-linearly on the gas temperature. Converting the ion density to neutral density 
is somewhat uncertain due to the unknown depletion factor and the fraction of molecular 
hydrogen in the CRL forming region. Based on the examination of the models for W48A 
and our experience in modeling other CRL data set (Roshi et al. 2005; Roshi et al. 2006), 
it is expected that the values for temperature in Table [1] can vary by a factor of 2 and the 
values for neutral density and L\\ can vary by a factor of 4. The model parameters are 
derived using a simplified geometry of plane parallel PDR slabs. A comparison of the PDR 
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properties obtained from such models toward W48A with those obtained from models which 
take into account the photo and chemical processes in the PDR shows that detailed modeling 
gives values within the factors quoted above (Jeyakumar et al. 2007 in preparation). 



3. Alfuen speed in molecular clouds and non-thermal width of carbon 

recombination lines 

Molecular clouds are largely neutral with a small admixture of partially ionized heavy 
elements with their free electrons (ionization fraction 10 -5 ). Any wave motion of the mag- 
netic fields in molecular clouds are strongly coupled to the ions. For those waves with A 
~ ^min, the time scale of momentum transfer between ions and neutrals is smaller than the 
magnetic perturbation time scale and hence neutrals are also coupled to such waves (Arons 
& Max 1975). For typical magnetic fields of 1 mG near 'dense' regions in star forming clouds 
and neutral densities of 10 6 cm -3 (eg. Johnston, Migenes & Norris 1989), X min ~ 3 x 10~ 7 
pc, which is much smaller than the size of these dense regions (fraction of a parsec). Thus 
magnetic waves with A <; 3 x 10 -7 pc can exist in dense regions and the characteristic speed 
(Alfuen speed) of these waves is determined by the total density (ie ion + neutral density) 
of the molecular cloud (Arons & Max 1975). 

Crutcher (1999) compiled the available sensitive Zeeman measurements of magnetic field 
strengths in molecular clouds and their neutral densities. Using this data, I estimate the 
Alfyen speed Va (units of cm s -1 ) in these clouds; 

Va= /A 2Blos , (1) 

where n# is the hydrogen atom density in units of cm -3 , \i = 1.4 is the effective mass 
of an H+He gas with cosmic abundance and mg is the mass of the hydrogen atom in gm. 
The measured line-of-sight magnetic field, Bi os , in units of G, is multiplied by 2 to convert 
it into total field strength (Crutcher 1999). Fig. [2] shows the estimated Alft;en speed in 
the molecular cloud sample taken from Crutcher (1999). The estimated Alfwen speeds are 
confined to a narrow range between 0.7 and 4 km s _1 with a median value of 1.6 km s _1 . 
This "constancy" of Alfven speed was noted earlier and has been understood from models 
of magnetic confinement of molecular clouds (eg. Moucschovias & Psaltis 1995; Basu 2000). 

We now compare the non-thermal velocity dispersion of carbon lines, a^T, observed 
from PDRs with Alfuen speed in molecular clouds. The non-thermal width is obtained by 
removing the thermal contribution, estimated using the gas temperature T PDR (see Tabled]), 
from the observed CRL width. The non-thermal line widths are plotted in Fig. [2] for the 
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corresponding PDR densities inferred from carbon line modeling. The non-thermal widths 
have values between 1.6 and 6 km s _1 with a median value of 2.9 km s _1 . Fig. [2] shows that 
the non-thermal line widths are almost "constant" over three orders of magnitude in density. 
The "constancy" of the non-thermal line widths and their magnitudes are similar (less than 
a factor of 2) to those inferred for Alfuen waves in molecular clouds. 

The gas phase carbon is ionized in PDR and hence is strongly coupled to the mag- 
netic field in these regions. If we consider similar physical parameters in PDRs as in dense 
molecular regions (n H ~ 10 6 cm -3 ; B ~ 1 mG) then X min of ~ 3 x 1CT 7 pc is at least two 
orders of magnitude smaller than the typical line-of-sight thickness of the PDR (ie Ly). Thus 
Alfuen waves with A <^ X m in exist in the PDR and here I consider their contribution to the 
non-thermal width of the observed carbon lines. The amplitude of the velocity of carbon 
ions, 8v, due to these waves is related to the magnetic perturbation amplitude 5B through 
the equation (eg. Arons & Max 1975) 

-■■ - SB (2) 

The observed non-thermal velocity dispersion is approximately given by 5v. It is usually 
assumed that 5B ~ B, in that case, the right hand side of Eq. (2) becomes identical to 
that of Eq. ([T]). Based on these considerations and the characteristics exhibited by the non- 
thermal velocity dispersion and Alfyen speed in molecular clouds (see above) we assume that 
the velocity dispersion of carbon lines is an estimate of the Alt-yen speed in the PDR. 

The observed velocity dispersion is usually scaled by v3 to convert it to 3D velocity 
dispersion. This scaling assumes random magnetic field orientation along line-of-sight and 
random polarization of Alfuen waves. Carbon lines are observed from PDRs near UCHs 
where shocks are present. In such shocked regions only the tangential component of the 
magnetic field is amplified and the Alfven speed associated with this component is scaled by 
the square root of the density compression ratio (McKee & Zweibel 1995). Hence the scaling 
factor needed to convert the observed velocity dispersion to a 3D dispersion is uncertain. 
Direct observation of the Zeeman effect of CRLs may help in determining this factor (see 
Section EJ). Here I note a systematically high value (a factor of 1.8) for the CRL velocity 
dispersion compared to the Alfwen speed in molecular clouds, which may be an indication of 
higher Alfuen speeds in PDR shocks. 



4. Magnetic field in Photo-dissociation regions 



The magnetic field strength is obtained using Eq flTJ by substituting the estimated 
non-thermal velocity dispersion of CRL for the Alfwen speed and using the neutral density 
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obtained from modeling the CRL emission (eg. Roshi et al. 2005b). Field strength values 
thus obtained are tabulated in Table. [TJ These values represent the total magnetic field 
strength in the PDR. Based on the expected range of the derived physical properties of the 
PDR (see Section [2]), the estimated uncertainty in B is typically a factor of 2.5. 

In Fig. [3j we compare the estimated magnetic field with those measured toward molec- 
ular clouds. The ordinate of the plot is the number density of H 2 molecules. Here, the 
estimated field strength values in the PDR are compared with those measured in molecular 
clouds with similar density. Such a comparison is possible since earlier observations toward 
molecular clouds show that the magnetic field scales with density (uh 2 oc p a47 ; Crutcher 
1999). To produce Fig. [3j the neutral density of the PDR given in Table. [TJis divided by 2 
to convert it into number density of H 2 . As seen in the figure the estimated magnetic field 
strengths in the PDR compare well within errors with those measured toward molecular 
clouds with similar density. 

Magnetic field measurements using Zeeman effect of CRLs near 1.4 GHz were attempted 
toward a few Hll regions (Silverglate 1984) of which W48 and S88B are of interest here. The 
1.4 GHz CRL emission toward W48 does not originate from the PDR associated with the 
UCH (Roshi et al. 2005a) and hence a comparison of the upper limit on the field strength 
obtained by Silverglate (1984) with the estimated value here is not meaningful. Toward 
S88B, the upper limit for the magnetic field strength obtained is consistent with estimated 
values given in Table [TJ Magnetic field measurements using OH or H I Zeeman effects are 
available in literature toward a few UCHs listed in Table [TJ I compare the field strength 
given in Table [TJ with the results of these observations. Note that OH and H I lines may 
originate from different spatial locations compared to the CRL forming region. From OH 
Zeeman observations toward S88B a field strength in the range 0.1 - 0.3 mG was obtained 
by Sarma et al. (2006 in preparation) consistent with our estimate. Van der Werf & Goss 
(1990) using Hi Zeeman observations measured a peak magnetic field of 0.1 mG in the —45 
km s~ 1 component observed toward W3. The CRL LSR velocity is comparable with this 
velocity and the field strength is consistent with the estimated value. Brogan & Troland 
(2001) measured a maximum field strength of 0.3 mG toward W49A. They used Zeeman 
effect of H I to measure the field strength with an angular resolution of 25". The measured 
value is about twenty times smaller than the estimated value; possibly because the two 
tracers (Hi and CRL) do not probe the same region in this case. 
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5. Discussion 

As mentioned in Section [21 CRL emission is dominated by stimulated emission for the 
transitions listed in Table [TJ Because of the dominance of stimulated emission, carbon line 
is detected only from the near side of the UCH. The non-detection from far side of UCH 
means that the width of CRL is not contributed to by the expansion of the UCH (if the UCH 
were expanding). Thus the line width has contribution only from thermal and non-thermal 
motions. The carbon line width may also have contributions from large scale motions such as 
outflows. However, the similarity of line widths observed from different sources (see Tabled]) 
may indicate that this contribution is small. Thus the field strength estimated in Section [4] 
may not be affected by such large scale motions. However, presence of any non-magnetic 
turbulence would affect the estimation of the field strength and hence the values obtained 
should be considered as upper limits. 

The magnetic origin of CRL width can be confirmed by measuring the field strength 
using Zeeman effect of carbon lines from the PDRs and comparing it with the values esti- 
mated in Section |H If confirmed then multi-frequency CRL observations form another tool 
to deduce the strength of magnetic fields near star forming regions. 

I acknowledge many useful discussions with K. Subramanian, S. Sridhar, A. A. Desh- 
pande, D. Bhattacharya and helpful comments from an anonymous referee. 
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Fig. 1. — Examples of carbon recombination line spectra toward UCHs G43.24— 0.05 and 
G70.33+1.59 observed near 9 GHz (Roshi et al. 2005a). The Gaussian profiles plotted in 
dotted line correspond to an assumed PDR gas temperature of 1000 K. The large difference 
between the observed line profiles and the Gaussian curves demonstrates the dominance of 
non-thermal motions in the PDR. The LSR velocity is with respect to the C89a (9.1779 
GHz) transition. 
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Fig. 2. — Alfuen speed in molecular clouds (crosses) and non-thermal velocity dispersion 
of carbon lines (triangles) observed toward PDRs near UCHs are plotted against density in 
these regions. Alfuen speeds are estimated from the molecular line data compiled by Crutcher 
(1999). Errors in the non-thermal velocity dispersion are derived from the uncertainties in 
the measurements and estimation of the model parameters. The expected range of the 
estimated densities for the PDR is shown by the horizontal bar. The plot shows that the 
non-thermal component of the CRL width is "constant" for almost all PDRs and has a 
median value similar to that of Alfuen speed in molecular cloud. 
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Fig. 3. — Log- log plot of the magnetic field strength vs neutral density. The plot compares 
the magnetic fields measured in molecular clouds (crosses; data taken from Crutcher 1999) 
with those estimated using recombination line data (triangles). Comparison is done for 
similar H2 densities in the molecular cloud and the PDR. Such a comparison is possible 
since earlier observations toward molecular clouds show that the magnetic field scales with 
density (Crutcher 1999). The plot shows that the field strengths estimated using CRL data 
(see Section. Hj) are comparable with those measured in molecular clouds with densities 
similar to PDR densities. The expected range of density and estimated magnetic field are 
shown by horizontal and vertical error bars respectively. 
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Table 1: Carbon recombination line data. 



Source 


Transition 


Ay 


TpDR 


n c + 


n H 


L|| 


B 


Ref 


name 




(km s^ 1 ) 


(K) 


(cm -3 ) 


(x 10 6 cm" 3 ) 


(x 10~ 3 pc) 


(mG) 




W48A 


C76a 


4.5(0.8) 


500 


2500 


8.5 


0.3 


2.9 


1 


W49G a 


C91a 


9.5(0.8) 


500 


3000 


10.0 


0.1 


6.8 


2 


W49J a 


C91a 


15.1(1.5) 


500 


2800 


9.0 


0.4 


10.4 


2 


S88B-east 


C92a 


4.4(0.3) 


600 


80 


0.4 


160.0 


0.6 


3 


S88B-west 


C92a 


5.6(0.4) 


400 


40 


0.2 


220.0 


0.6 


3 


GGD12-15 


C92a 


4.7(0.7) 


330 




0.02 


130.0 


0.1 


4 


W3A 


C168a 


5.5(0.9) 


100 




0.1 


50 


0.4 


5 


G32.80+0.19 b 


C89a 


7.5(1.3) 


500 


1000 


3.4 


1.1 


3.1 


6 


G37.87-0.40 b 


C89a 


14.0(1.5) 


500 


120 


0.4 


29.0 


2.0 


6 


G43.24-0.05 b 


C89a 


6.2(0.9) 


500 


300 


1.0 


9.0 


1.4 


6 


G45.12+0.13 b 


C89a 


12.9(1.9) 


500 


40 


0.1 


124.0 


0.9 


6 


G45.45+0.0 b 


C89a 


9.6(2.9) 


500 


30 


0.1 


150.0 


0.7 


6 


G70.29+1.60 b 


C89a 


8.3(0.5) 


500 


350 


1.2 


8.6 


2.1 


6 


G70.33+1.59 b 


C89a 


4.2(0.4) 


500 


170 


0.6 


25.0 


0.7 


6 



a Modeling of carbon line emission for these sources has provided upper limits for n c +, nn and lower limits 
for Lii . 

h The data for these sources are taken from the Arecibo survey. Modeling of carbon line emission for these 
sources has provided lower limits for nn and upper limits for L». 

References: (l)Roshi et al. (2005a); (2) Roshi et al. (2006); (3) Garay ct al. (1998); (4) Gomez et al. (1998); 
(5) Kantharia et al. (1998); (6) Roshi ct al. (2005b). 



